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Abstract 

Based on static and dynamical density functional theory, a phase-field-crystal model is derived 



which involves both the translational density and the orientational degree of ordering as well as 
a local director field. The model exhibits stable isotropic, nematic, smectic A, columnar, plastic 

00 

CN crystalline and orientationally ordered crystalline phases. As far as the dynamics is concerned, 

the translational density is a conserved order parameter while the orientational ordering is non- 
§ conserved. The derived phase-field-crystal model can serve for efficient numerical investigations of 

various nonequilibrium situations in liquid crystals. 
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I. INTRODUCTION 



Within the phase-field-crystal (PFC) model [IHS] the crystalline density field is described 
basically in terms of a single Fourier mode, i.e. as a sinusoidal density wave. The PFC 
model can be understood as a modified Landau expansion of the full inhomogeneous one- 
particle density field of the solid. It has been applied for large-scale numerical investigations 
of dynamics in the solid state. Characteristic examples include calculations of a variety of 
different quantities: the structure and free energy of the fiuid crystal interface [HIS], crystal 
growth dynamics into a supercooled liquid [H], the structure [7] and dynamics [S] of grain 
boundaries, and the Asaro-Tiller-Grinfeld instability [SHIS]- A solid particle just enters as 
a "blob" , i.e. a weak density modulation, and the dynamics is diffusive on long time-scales, 
i.e. the density itself is a conserved order parameter. 

Recently, the PFC model was derived from dynamical density functional theory (DDFT) 
[T3HT6] . Static density functional theory provides a microscopic framework to describe crys- 
tallization in equilibrium piTH^ and a Landau expansion in terms of density modulations 
[2TH21] can be used to derive the corresponding approximative free energy for the PFC. 
Density functional theory was generalized towards nonequilibrium dynamics for Brownian 
systems [251 - I27] and the resulting dynamical density functional theory can be used to derive 
the dynamics of the PFC [16]. First of all, this derivation should apply to colloidal disper- 
sions whose short-time dynamics is clearly diffusive. But also molecular systems governed 
by Newtonian dynamics for short-times behave diffusive on longer time-scales and therefore 
the derivation might have relevance for atomic systems as well. 

The PFC model has been generalized to mixtures by including more than a single density 
field |23]. However, it has never been applied to liquid crystals which are made by particles 
with orientational degrees of freedom [28j. Under appropriate thermodynamic conditions, 
these particles occur in liquid-crystalline phases including nematic, smectic A, columnar, 
and plastic crystalline phases [2^ [5U] . While the DDFT approach was recently generalized 
towards orientational dynamics for Brownian rods both in three [31] and two [32] dimensions, 
the link towards the PFC model has not yet been elaborated for orientational degrees of 
freedom. 

In this paper, we close this gap and propose a PFC model for liquid crystalline phases. 
One motivation here is to propose a minimal model, i.e. the simplest nontrivial model for 
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dynamics of liquid crystals. We derive this model from dynamical density functional the- 
ory. Depending on the model parameters, the resulting model does accomodate isotropic, 
nematic, smectic A, columnar, plastic-crystalline phase and an orientational ordered crys- 
tal. It can therefore be used to describe the statics and nonequilibrium dynamics in various 
situations where these phases are relevant. This may stimulate further numerical investi- 
gations. The model is basically formulated in terms of two density fields, a translational 
and an orientational one, plus a local director field. While the translational density field is 
conserved, the orientational one and the director field are nonconserved and relax quicker. 

The paper is organized as follows: in section II, we derive the PFC model from fluid- 
based density functional theory by expanding the orientational dependence of the density 
field up to the first nontrivial order and performing a gradient expansion in the translational 
coordinate. Then, in section III, we discuss the parameter space for which stability of the 
different phases is obtained. The dynamical equations are derived from dynamical density 
functional theory in section IV. We then discuss possible extensions of the model to more 
complicated situations and give final conclusions in section V. 

II. DERIVATION OF THE PHASE-FIELD-CRYSTAL MODEL FOR LIQUID 
CRYSTALS: STATICS 

We start our derivation from microscopic static density functional theory for liquid crys- 
tals. We consider particles with orientational degrees of freedom described by a set of 
unit vectors {ui;i = 1, A^} and center-of-mass positions {Ri;i = 1, A^}. Though most 
of the considerations can be done in three-dimensional space, we restrict ourselves in the 
following to two spatial dimensions, where Ri G and Ui{(f)) = (cos 0, sin 0) (0 G [0,27r[). 
The system has a total area A and is kept at finite temperature T. 

A pair interaction potential V{Ri — R2,ui, U2) between two particles 1 and 2 is assumed. 
We henceforth consider apolar particles implying the following symmetries 

V{f,Ui,U2) = V{-f,Ui,U2) = V{f, -Ul,U2) = V{r,Ui, -U2) = V{f,U2,Ul) (1) 

Examples for V{Ri — R2,ui,U2) comprise: i) excluded volume interactions as dictated by 
hard spherocylinders [301 ES] or hard ellipsoids |3l] which are used for sterically-stabilized 
colloids, ii) Yukawa segment models |3SHSZ] used for charged colloidal rods, iii) Gay-Berne 
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potentials j38] - li0] used for molecular liquid crystals. 

The inhomogeneous one-particle density p{R, u) provides the joint probability density to 
find particles at center-of-mass-position R with orientation it. 

p{R, u) = ( ^ 6iR - R,)5{<P - (2) 



i=l 

where for an observable A 

f.27r i'2-K 



{A) = ]- I SRi... / cP'Rn I d(j)i... I d(j)N Aexp 



VjRi - Rj,Ui,Uj 



(3) 

A Jo 

is the normalized canonical average, fcs denoting Boltzmann's constant, and the classical 
canonical partition function Z ensures the normalization < 1 >= 1. Clearly, for apolar 
particles, p{R, u) = p{R, —u). 

Classical density functional theory of inhomogeneous fiuids now provides the existence of 
an excess free energy density functional such that the functional 

p /'27r 

n(T,A,/i,[p(i?,M)])= J-,rf(T,A,[p(i?,M)])+J-,,.,(T,A[p(i?,ti)])- / d^R / d<Ppp{R,u) 

J A Jo 

(4) 

is minimal for the equilibrium density field for a given chemical potential p, temperature T 
and area A. The ideal rotator gas functional J^id is known exactly: 

T,diT,A,[p{R,u)]) = kBT / d^R / c?0p(i?,M(0))[ln(AV(^,M(0)))-l], (5) 

J A Jo 

where A denotes the (irrelevant) thermal wavelength. The excess free energy functional 
J^exc{T,A, [p{R,u)]), on the other hand, incorporates all correlations and is not known in 
general. In the low density limit, a second virial approximation (Onsager functional) is 
getting asymptotically exact [29] 

2n r^TT 



J^,,,{T,A,[p{R,u)])^l [ d'R, [ d^R2 r d<P, f 

^ J A J A Jo Jo 



X I exp I J^^Rli:^}^ , _ 1 , p^R,^uMR,,u,). (6) 

More generally, the Ramakrishnan-Yussouff theory of freezing [T7] can be applied to get the 
following perturbative approximation for J^exc 



i02 



j-e,e(T,A[p(^,^)])~-^ / d'^Ri I d'R2 r d<p^ r 

o A J t/ •-'0 

X c^'^\Rx - %MM){p{R\M) - p){p{R2M) - p) (7) 
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which can be viewed as a truncated density expansion in the density difference p{Ri, ui) — p 
around a mean density p with the Icernel representing the direct correlation function of the 
reference fluid at temperature T and density p. 

Another expression which works complementary at high density for very soft interactions 
[21] is a mean-field approximation 

Fexc{T,A,[p{R,u)]) ^ - / (fRi / £R2 I d(f)i / d(f)2V{Ri - R2,Ul,U2)p{Rl,Ui)p{R2,U2) 

^ J A J A Jo Jo 

(8) 

More accurate forms for J-'exc have been proposed for hard particles using weighted-density- 
approximations im 132] or fundamental-measure theory [20] . 

In the following we shall adopt the Ramakrishnan-Yussouff theory and approximate fur- 
ther by only considering weak anisotropics in the orientations. The leading expression in 
the density parametrization is then 

piR, u)=p + pV^i(f) + p^2(r) ({u ■ Mo(f))2 - + - (9) 

Here, the real- valued dimensionless orient at ionally averaged density is ipii"^) which is iden- 
tical to the original treatment of the PFC model [H [2]. The dimensionless field ^'2(^5 on 
the other hand, measures the local degree of orientational order. For apolar particles, the 
leading anisotropic contribution is the third term on the right-hand-side of Equ. ([9]). Finally, 
the field UQ{f) defines the local director of the orientational field. [13] 
We now derive the static free energy functional. With 

X = iJi + ^2P2{U ■ Uo) (10) 

where P2{y) = — \i the ideal rotator gas part reads as 

7id = ksTp / d^R / d0 (1 + x) [In {K^p{l + x)) - l] 

J A Jo 

= Fo + kBTp J^^'r/J d0 Qx^ - ^x' + ^x' + 0{x')^ (11) 
where Fq = 2TTAkBTp(\n{A'^p) — 1) and irrelevant terms linear in x on the right hand side 



of Equ. (11) were absorbed in a scaled chemical potential. Inserting (10) and performing 



the angular average, we obtain 



(12) 
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The correlational part within the Ramakrishnan-Yussouff approximation is 



2lT 



2-IT 



1 / d02hAl(i?l)+V^2(i?l)P2U(0l)-Mo(i?l) 



A J A Jo ~'0 

^l(i?2) + V'2(i?2)P2 fw(02) ■ Wo(i?2))) C^^H^l " ^2, 01, ^2) 



(13) 



We now decompose 



00 00 



c(^,0i,02)= J2 cw(i?)e'''"^^e2*™'^^ (14) 

m=—oo m'=— 00 

and consider only the leading terms where m,m' G { — 1,0,1}. The relevant expansion 
coefficients are 

r2n r2n 



(R) 



(27r)2 



d(f)' e 



I -2im4>-2im'(l)'j2) 



(15) 



'0 JO 

By symmetry, it can be shown that Coo(-R), c_ii(i?) and Ci_i(-R) only depend on \R\. There- 



fore 

•F txc 



ksTp 



d?Ri I d^R2^7c^ 



cooilRi - R2\)MRl)MR2) 



+ -MRi)MR2) [co-iiRi - i?2)e-'*'^°(^^) + Coi(i?i - i?2)e'^'^«(^^) 
+ \mR2)MRi) {c^ioiRi - 4)e-2^'^°(^^) + cio(^i - ^2)e^^^«(^^) 
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MRi)MR2) {c-i-i(i?i - R,)e-^^MR.)-2^MR2) + ^^^^q^^ _ 



2i(t>o{Ri)+2i<Po{R2) 



2iMRl)+2i<t>o{R2) 



(16) 



Now a gradient expansion is performed [23] up to fourth order in the ipiipi term of Eqn. (16) 
and up to second order in the 1^11^2 and V'2'0'2 terms. We assume that the highest gradient 
term ensures stability. Thereby one obtains 



exc exc ^ exc ~ exc 



(17) 



with 



and 



77(1) 

exc o 2 - / 12 D 



^ \ 2 



A^f(/2)-S V^i(i?) +C A^i(i?) 



77(2) 

ex 



exc r, 2- / 12 



27r^p / d'R 
ksT J A 



Dijl{R) + E{(ViJ2{R)]\^'il^l{R) i^MR)^ 



(18) 



(19) 



and 

7r(3) 
•> exc 



A 



2Ti'-p d'RF VMR) ■ v^2(i?) + 2^/j2iR)iMR) ■ WMR) ■ (20) 



In detail, in (20), {uo{R) ■ V)^ := Y.- ■^Uoi{R)uQj{R)didj where Moi(-R) 



cos (j)o{R) 
sin0o(-R) 



In (18)-(20), the coefficients A, B,C, D, E and F are generalized moments of the direct 
correlation function. In general, they depend on the thermodynamic conditions (T, p). In 
detail [S] (for A = M^)^ 



oo 



A = -27rp / dR RcooiR) (21) 







oo 



B = Tip dR R^coo{R) (22) 



oo 



C = ~l dRR%o{R) (23) 







oo 



D = -!^ / dRRc.uiR) (24) 
E = Tip 







POO 

/ dRR^c.niR) (25) 
Jo 

poo 

F = -- dR R^ coi{R COS (j)R,R sin ^B)e'^'''"'. (26) 
8 Jo 

As a remark: F does not depend on 0r. For stability reasons, we henceforth assume 
C,E>0. 

Let us now discuss the static free energy functional. In the limit of no orientational 
order, 1^2 = 0, one recovers the phase-field crystal model of Elder and coworkers PQ The 
expansion up to fourth order is formally similar to a Landau expansion of the smectic A- 
isotropic phase transition if ipi represents the smectic order parameter |15] . In the opposite 
case of constant ipi and constant 1^2, Frank's elastic energy with a nonvanishing splay and 



vanishing bend modulus is recovered in the term ~ (V0o(-R))^ iii (19). In fact, in two spatial 
dimensions there are only two Frank elastic constants since the twist modulus vanishes. If ipi 
is constant and both ip2 and Uq are space dependent, we obtain the Landau-de Gennes free 
energy |16] used for inhomogeneous uniaxial nematics. In two dimensions, the Landau-de 
Gennes free energy has only one gradient coefficient |47j which is again the parameter E. 
For space-dependent ipi and uq, but constant ^'2, the free energy derived by Pleiner and 
Brand [H] is recovered. Finally, for space-dependent ipi and 1^2, we obtain the coupling 
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liquid crystalline phase 




"02 


Uq 


isotropic 








irrelevant 


nematic 


constant 




constant 


plastic crystalline 


oscillatory 





irrelevant 


orient ationally ordered crystalline 


oscillatory 




constant or oscillatory 


smectic A 


planar oscillatory 




constant, oscillatory V^i 


columnar 


planar oscillatory 




constant, oscillatory _L Vipi 



TABLE I: Characteristic values for the number density -01, the nematic order parameter -02 and 
the director field uq for six different liquid crystalline phases, namely isotropic, nematic, plastic 
crystalline, orientationally ordered crystalline, smectic A and columnar. 



terms in J^exc proposed by Brand and Pleiner |49J. However, the full free energy functional 



(17) with fourth-order gradients in ipi and the appropriate couplings to ip2 and uq is new 



and constitutes the basic static result of this paper. 



III. EQUILIBRIUM BULK PHASE DIAGRAM 

By minimizing the free energy functional, for given thermodynamic parameters T and 
p, the equilibrium phase diagram is gained. In the special case of tp2 = 0, the PFC phase 
diagram of Elder and coworkers [2] is obtained. By scaling out a length scale, there are 
only two remaining parameters for which a fluid, a triangular phase and an (unphysical) 
stripe phase is stable (see Fig. 4 in [2j). For D < — |, a nonzero stable value for iIj2 
occurs. Combined with the PFC phase diagram, the possibility of a nematic phase and 
an orientationally ordered crystal emerges. The stripe phase at ^/'2 7^ becomes either a 
smectic A or columnar phase depending on the sign of the parameter F. In fact, it was 
already shown in Ref. [50j that an Onsager-functional yields a smectic A phase. 

All the possible liquid crystalline phases are summarized in Table [T] together with their 
characterizing values for the number density ipi, the nematic order parameter ip2 and the 
director field uq. For D > — |, a plastic crystal and the ordinary isotropic phase can be 
stable. The full numerical calculation of the equilibrium phase diagram as a function of the 
parameters A, B, C, E, and F is planned in a future study. 
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IV. DERIVATION OF THE PHASE-FIELD-CRYSTAL MODEL FOR LIQUID 
CRYSTALS: DYNAMICS 



A. Dynamical density functional theory 

In two spatial dimensions, the dynamical density functional theory for Brownian systems 
is a deterministic equation for the time dependent one-particle density field u, t) 



dt \ keT 5p{R,u,t) ^0 1 d<P5p{R,u,t) 



(27) 



Here, is the diagonal translational short-time diffusion tensor which we assume to be 
isotropic in the following, Dt = diag{DT, Dt), and is the rotational diffusion constant. 
Furthermore, 7 = 7id + 7exc is the total free energy functional. If the density parametriza- 
tion ([9]) is used, this becomes a functional 7\il)i{R),il)2{R)-,4'o{R)\ of the three scalar fields 
■ipi{R), 'iIj2{R), 4'o{R)- Now the chain rule of functional differentiation yields: 

57 1 57 ,4 57 ^ rm^^ ■ ^ sin(2(0 - 0o(i?))) 

- + 7^-P2(cOS(0 - 0o(-R))) 



5p{R,(f)) 27rp5^^(i?) TipSip^iR) t^P5(Pq{R) 

(28) 



By inserting this into Eqn. (27), coupled equations of motion can be obtained. 



B. Derivation of the dynamics (PFCl model) 

First we describe the dynamics for the most case which is called PFCl model in Ref. 
16] . The PFCl model avoids two further approximations, namely the expansion of the 



logarithm (11) and a constant mobility assumption. By inserting the chain rule (28) into 



the dynamical density functional theory (27), one obtains dynamical equations for the three 
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scalar fields il)i{R,t),il)2{R,t) and (f>o{R,t) as follows: 



ksTpi,, = kBTpDTA^, + — (Iv f (1 + ^i) 



kBTp^2 = ksTp ( Dt A^2 - 4Dr^2 - 4DtV'2 (v0o)' 



5^2 



(29) 



+ - V 4(l + #,) V^-^^ +yV— 



+ (1 + V'l) 



-lb- '^-^ ^ 



^(V0o) -4V 



6o V's 



V0O 



D 



R 



16(1 + ^,)- 
(30) 



kBTp^24>o = ksTp (DtV^2A0o + • V0o)) 

4(1 + V'i 



+ 



+V (1 + V^i) 4 



^^2 



^V0o + V 



500 ^1 



TT 



4(l + ^i)^- 
d(po 1p2 

(31) 



The right-hand-side of Eqn. (29 ) clearly shows that the time-derivative ipi is proportional 



to a divergence of a generalized current. This implies that a generalized continuity equation 
holds such that the order parameter field ipi{R,t) is conserved. On the other hand, this 
is not true for the two remaining orientational order parameter fields ip2{R,t) and (j)Q{R,t) 
which are therefore non-conserved. 
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The functional derivatives are local and given by 



F ^ 

-Aij2 + fY, 9idj {^2UoiUoj) , (32) 



F 





2 



- -^^ipi + P 2^ uoiUojdidj^pi, (33) 

^ = -AE^lA<^o + F^P2 E (34) 
d(Po O(po 

where ^!fo^ ^ ^ 20o cos 20o\ _ ^3^^ 
^<^o \^ cos 200 sin 200/.. 

Combining these equations yields explicit deterministic and coupled equations of motion 
for the three order parameter fields ipi{R, t), '02(-R, t) and 0o(-R, t) which can be implemented 
for a numerical solution. 



C. Derivation of the phase-field-crystal model with constant mobility (PFC2 
model) 



In the constant mobility approximation, the prefactor in front of the density functional 
derivatives on the right-hand-side of Eqn. (27) is replaced by the constant -j^. Then the 



dynamical density functional equations simplify to 



dp{R, u, t) 
dt 



DtA + D 



6T 



kBTSp{R,u,t) 



(36) 



In this case, the equations of motion for the three scalar fields ipi{R, t),ip2{R, t) and 0o(-R, t) 
read as 



kBT-npipi = \dtA^ 
2 oipi 



kBT7rp-ip2 = Dt 



4A— - 16 V0o^^ - M ^ A;r^ ■ ^'^o - 4— 



(37) 
6T 



6lp2 

(38) 
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and finally 



8 V 



5^ 

5^2 



V0O + 4— A0O + A 



V0C 



-4D 



1 (5J^ 



^2 500 



(39) 



The ordinary phase-field crystal model is obtained by a subsequent expansion of the ideal 



rotator term (11) up to fourth order. Following Ref. |16j, the resulting dynamics is called 
PFC2 model. In this case, the density functional derivatives are again local and given by 



1 



1 



71 [2 + 2^1 -^f 



2 ,n ?/'l-0, 

3^^ 4 



p 



(40) 



vr 



4 4^^^' + t^'"^^ + 64 



+ An' ( Z}^2 - ^A^2 + 4EV^2 (^V0o 



57" 



fcsTp 500 



-4E^/'2A0o + F^2 



duoiUoj 



(41) 
(42) 



The advantage of these equations is that they reduce to the dynamics of the traditional 
phase-field-crystal model in the pure translational case. For a rough numerical exploration, 
the PFC2 model should give the same qualitative answer as the PFCl model. For spherical 



particles this was shown in Ref. [16]. The dynamical equations (37)-(39) represent the main 
result of this paper. 



V. CONCLUSIONS 



In conclusion, we derived from static and dynamical density functional theory phase-field- 
crystal equations which govern the diffuse nonequilibrium dynamics for liquid crystalline 
phases. The approximations involved a two-fold: first the density functional is approximated 
by a truncated functional Taylor expansion similar in spirit to the Ramakrishnan-Yussouff 
theory. Then a generalized gradient expansion in the order parameters is performed which 
leads to a local density functional. In addition to the traditional scalar phase-field variable 
■01, a local scalar nematic order parameter -02 and a local nematic director field 0o was 
introduced and coupled to the phase-field variable ipi. If the additional variables are zero, 
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the phase-field-crystal model of Elder and coworkers [H |2] is recovered. If, on the other hand, 
tpi is set to zero we recover the Landau-de Gennes free energy for uniaxial nematics extended 
by Pleiner and Brand jlHl US]. The proposed phase-field-crystal model for liquid crystals 
allows for a wealth of stable liquid crystalline phases including isotropic, nematic, smectic A, 
columnar, plastic crystalline and orient at ionally ordered crystals. How the stability of these 
phases depends in detail on the model parameters still needs to be explored numerically. 
The new coupled phase-field-crystal equations can be used to simulate the nonequilibrium 
dynamics of liquid crystals. Possible problems are dynamics of topological defects in the 
nematic phase [ST] and the formation of metastable phases at a growing interface [52]. As 
the dynamics in nematic states can be obtained by using other approaches like the one 
in Ref. [19], the present model may be applicable in particular to smectic films and to 
two-dimensional crystalline phases. 

In the present paper, the derivation of the phase-field-crystal model was performed in two 
spatial dimensions. Though more tedious there is no principle problem in doing the same 
analysis in three spatial dimensions with the use of spherical harmonics for the orientational 
degrees of freedom. Moreover the present derivation can in principle be done to higher 
order in the orientational degrees of freedom. The translational degrees of freedom can be 
anisotropic for the dynamical mobility matrix [53j. 
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